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We investigate the entanglement properties of an ensemble of atoms interacting with a single 
bosonic field mode via the Dicke (superradiance) Hamiltonian. The model exhibits a quantum 
phase transition and a well-understood thermodynamic limit, allowing the identification of both 
quantum and semi-classical many-body features in the behaviour of the entanglement. We consider 
the entanglement between the atoms and the field, an investigation initiated in Holl .In the ther- 
modynamic limit, we give exact results for all entanglement partitions and observe a logarithmic 
divergence of the atom- field entanglement, and discontinuities in the average linear entropy. 



I. INTRODUCTION 

Understanding entanglement - the quantum correla- 
tions impossible to mimic with local classical theories - is 
a fundamental goal of quantum information science. Sim- 
ilarly, understanding complex modes of behaviour, such 
as quantum phase transitions Q and quantum chaos 0, 
has become an important part of quantum many-body 
theory. Since large correlations and collective behaviour 
are an intrinsic part of critical systems, concepts and 
formalisms used to describe entanglement are now being 
employed to reveal the truly quantum nature of certain 
aspects of criticality. 

Investigations into the entanglement between interact- 
ing spin-1/2 systems on a one dimensional chain have 
revealed so-called 'critical entanglement', in which an 
entanglement measure of the ground state exhibits uni- 
versality, or scaling behaviour, around the critical point 
In particular, for the infinite XY spin chains, (and 
their Ising variants) it has been shown that entanglement 
between nearest and next-nearest neighbours reaches a 
maximum near, but not at, the critical point 0,0. Fur- 
thermore, Osterloh et al. Q have observed scaling be- 
haviour of the entanglement, showing that the derivative 
of the concurrence diverges logarithmically near the crit- 
ical point. They also found a logarithmic divergence of 
the derivative as a function of system size. 

Latorre, Vidal, and co-workers [f|, |i| took a differ- 
ent approach and investigated the entanglement, via the 
von Neumann entropy, between a block of L spins and 
the rest of the chain in XY and Heisenberg spin-chains. 
They found a logarithmic scaling of the entropy with L; 
this time with a pre-factor corresponding to the 'central 
charge' of a 1 + 1 continuum quantum field theory of the 
same universality class. In effect, they found the same 
area law associated with the geometric entropy first stud- 
ied by Srednicki Q . In an effort to understand the nature 
of the scaling of entanglement, Orus et al. @ illustrate 
that the scaling of the entanglement at the critical point 
determines whether or not one could efficiently simulate 
the quantum system at this point on a classical computer. 

Going beyond 1-dimcnsional spin chains, the authors 
of [j| studied a highly connected simplex, where each 



spin interacts equally with all other spins, and the lattice 
spacing no longer plays an important role. Importantly, 
because of the symmetry, they find a maximum in the 
pairwisc concurrence at the critical point, and determine 
scaling exponents for the behaviour of the concurrence 
with system size. 

In this paper, we continue the investigation, begun 
in [lfjj . of the entanglement properties of the single- 
mode Dicke Hamiltonian, which describes an ensemble 
of N two-level atoms coupled to a single-mode bosonic 
field. This model exhibits a 'supcrradiant' quantum 
phase transition (QPT) in which the ground state un- 
dergoes a dramatic change in character. We consider 
several aspects of the ground-state entanglement in this 
model and observe how they are affected by the QPT. 
We investigate entanglement between the atomic ensem- 
ble and the field mode via the von Neumann [Til IT3 | and 
linear entropies of this bipartite decomposition. We also 
calculate the average linear entropy of all the subsystems, 
which corresponds to a multipartite measure introduced 
by Meyer and Wallach [T^]. In the thermodynamic limit, 
the model is exactly soluble across the whole coupling 
range, and we give exact results for these quantifiers of 
the entanglement. For finite N we use perturbative and 
numerical methods. 

In summary, we find the atom-field entropy diverges 
at the phase transition alongside the traditional correla- 
tion length, with corresponding critical exponents, and 
may be fruitfully described by an effective 'entanglement 
temperature'. As has been discussed previously 
the QPT is foreshadowed at finite N by various 'precur- 
sors', and in particular, a transition from integrable to 
Quantum Chaotic behaviour near the critical point. This 
transition is characterised by a change in the energy level 
statistics, and can be correlated with the change in the 
phase-space of a classical Hamiltonian corresponding to 
the Dicke model. The phase transition in the quantum 
model maps to a supercritical pitchfork bifurcation in the 
classical model, and such bifurcations have recently been 
related to entanglement characteristics 0| . In ad- 
dition, Fujisaki et al. have shown that the appearance 
and strength of chaos can be linked to the production of 
entanglement |Tsj . Further work is required in clarifying 
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the relation between entanglement in quantum systems 
and chaos in the corresponding classical model. However, 
there is a conceptual connection between the divergence 
of trajectories in classical chaos and the derealization of 
the quantum ground state, which is, in general, indicative 
of entanglement. 

The model considered here is of wider interest still, 
given that the interaction of a charge or spin systems 
with a single bosonic mode is viewed as a mechanism for 
generation of entanglement in many different situations 
such as quantum cavity QED, quantum dots [T^. l20l |. 
and ion traps. In addition, many suggestions have been 
made to use the environment, or bosonic cavities, to share 
or mediate entanglement [13, I^H |23| ■ In particular, 
Reslen et al [24| have shown that there is a direct equiv- 
alence between the single mode Dickc Hamiltonian and 
the infinitely-coordinated XY model. 

This paper has the following structure. In section ITT1 
we reintroduce the Dicke Hamiltonian, and describe the 
quantum phase transition. In section ITTT1 we consider the 
atom-field entanglement by recalling the finite numeri- 
cal and exact thermodynamic limit results for the von- 
Neumann entropy considered in [Tol |. and, as mentioned, 
extending the discussion with a calculation of the linear 
entropy, participation ratio, and the average linear en- 
tropy. We omit discussion of the pairwise entanglement 
covered in [T^j. We conclude with discussions in section 

II. THE DICKE MODEL 

Generically, the Dicke Hamiltonian (DH) describes the 
dipole interaction between TV atoms and n bosonic field 
modes. Here we shall only consider the single mode case 
with n = 1. A standard approach to such quantum-optics 
Hamiltonians is to make the rotating wave approximation 
(RWA), rendering the model integrable. We do not make 
the RWA here, allowing the model to describe both weak 
and strong coupling regimes. 



There exists a conserved parity operator 

n = e < *t° to+J * +J '> > (2) 

which commutes with the Hamiltonian. For finite N, the 
ground state has positive parity. The DH undergoes a 
QPT at a critical value of the atom-field coupling A c = 
y/WLUQ /2 which breaks this symmetry. 

At finite TV, we perform numerical diagonalisations us- 
ing a basis \n) ® \j, m), where |n) are Fock states of the 
field, and \j,m) are the so-called Dicke states - eigen- 
states of J 2 and J z . We make use of the parity symmetry 
to simplify these numerics. 



B. Thermodynamic Limit 

The DH undergoes a QPT in the thermodynamic limit 
(j — > oo N — > oo, notation which we will use inter- 
changeably) at a critical coupling of A c = ^/ujujq/2. Be- 
low A c the system is in its normal phase in which the 
ground state is largely unexcited. Above A c , the super- 
radiant phase, the ground-state possesses a macroscopic 
excitation. 

As illustrated in Ref. [l4|, exact solutions may be ob- 
tained for both phases in the thermodynamic limit by 
employing a Holstein-Primakoff transformation of the an- 
gular momentum algebra. In this section, we briefly sum- 
marise this analysis, highlighting those features such as 
are required here. 

The Holstein-Primakoff mapping expresses the angular 
momentum in terms of a single boson mode, 

J+ = 6 f y/2j - tfb, J- = y/2j - tfbb, J z = tfb- j,(3) 

with [b, b'] = 1. In this representation, the DH trans- 
forms into a two mode bosonic problem. 

1. Normal Phase 



A. The Hamiltonian 

The single-mode Dicke Hamiltonian is 

N N . 

i=l i=l » " 

= lu q J z + uja^a -\ ^= (a 1 + a) (J + + J_) , (1) 

where J z = X)i!=i s zj J± = Sj=i s ± are collective an- 
gular momentum operators for a pseudo-spin of length 
j = N/2. These operators obey the usual angular mo- 
mentum commutation relations, [J Z ,J±] = i</± and 
[J+, J-) = 2J Z . The frequency w describes the atomic 
level splitting, u> is the field frequency, and A the atom- 
field coupling strength. 



The normal phase is found by simply taking j ^ oo 
in the bozonised hamiltonian, which produces a linear 
two mode hamiltonian. This, as described in [Toj . can be 
diagonaliscd with a Boglioubov transformation. 

To calculate the atom-field entanglement of the ground 
state, we require the reduced density matrix (RDM) of 
the atoms in the ground state. Summarising our steps in 
|10( , the ground-state wave function is a product of two 
Gaussians, ^(gi,^) = G + (qi)G-(q2) described by the 
co-ordinates corresponding to the bosonic operators of 
the diagonaliscd Hamiltonian. Inverting the Boglioubov 
coordinate rotations gives us the wave function in terms 
of the coordinates (x, y) corresponding to the physical 
field (x) and atom (y) modes. To obtain the RDM of the 
atomic system, we integrate over the x coordinate. We 
write the resulting RDM in terms of a rescaled y coor- 
dinate y — > y/n (writing c = cos7^ 1 -', and s = siii7 < - 1 ^, 
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Pa(y,y') 



7r(e_c 2 + e+s 2 ) / 
2e_e, +D 



x exp 



4k 2 (c_c 2 + e+s 2 ) 
D 



(if 



2« 2 (e_c 2 + e+s 2 ) 



2/2/ 



(4) 



where D = (e_ — e+) 2 c 2 s 2 , and the excitation energies in 
this normal phase are 

4 = 5 ( w o + ^ 2 ± \A w o - ^ 2 ) 2 + 16A 2 ^ ) • (5) 

We did not perform this rescaling in [ljj ■ As this rescal- 
ing is effected by a unitary transformation on the atomic 
system alone, it will not affect the atom-field entangle- 
ment. It does, however, aid in the interpretation of our 
results, as we show later. Note that the RDM for the 
field mode is the same as above, except with c and s 
interchanged. 

2. Super-Radiant Phase 

In the following section, we describe in more detail the 
calculations and properties of this phase which were not 
covered in ^(J- I n the superradiant (SR) phase (A > A c ), 
both atom and field degrees of freedom acquire macro- 
scopic mean-fields. We incorporate these mean-fields by 
displacing the two oscillator modes 



,t 



± V», &t = S T y/fl, 



(6) 



where a, (3 are of order j. That there are two choices 
of sign here is significant, as the two choices lead to two 
different Hamiltonians with degenerate solutions - an in- 
dication that the parity of the system has been broken 
in this phase. 

By inserting one of the above displacements into the 
Holstein-Primakoff bozonised Hamiltonian, and setting 
terms with overall powers of j in the denominator to 
zero, we obtain an exactly soluble Hamiltonian. Diago- 
nalization requires a specific choice for the displacements 



= 17 VK 1 - M), = 0'd - M), with n = A 2 /A 2 , 



UJ V 2 

and a rotation of the coordinates 

Qi = Xcos 7 (2) - Fsin 7 (2) , 
Q 2 = Xsin 7 (2) +ycos 7 (2) , 



(7) 



with angle of rotation given by 

tan(2 7 ^) = 4^. (8) 
The excitation energies of the SR phase are 



and e_ is real only for A > A c . The effective Hamiltoni- 
ans derived with either choice of sign in Eq. (jfj do not 
commute with the parity operator n, and thus we see 
that this symmetry is broken in the SR phase. 

As before, the ground state of the diagonalised Hamil- 
tonian is the product of two Gaussians in Q± and Qi- To 
obtain the wave function in terms of the original atomic 
and field co-ordinates, we must not only perform the ro- 
tation Qi, Q2 — » X, Y but also take into account the rela- 
tionship between the displaced and re-scaled coordinates 
X, Y and the original atom-field co-ordinates x, y, 



X = xt 




(10) 



where uj = t^(1 + m)- I 11 the displaced (X, Y) frame 
the wave functions have the same form as in the normal 
phase, but with different parameters and coefficients. In 
the original frame (x, y) the wave functions are again the 
same but displaced from the origin. 



3. Two Lobes 

The two possible displacements lead to two Hamilto- 
nians with ground state wave functions in the x, y repre- 
sentation displaced from the origin in different directions. 
For large but finite N, the ground state in the SR phase 
is a single two-lobed wave function, 



SR/ 




—{y±\—)s 

UJ V w 



2a, 



UJ V 



(11) 



where G± are the normalised Gaussians. This state has numerical results at finite -/V: since the displacements 
positive parity and can be used for comparison with our are of order y/N, for N 3> 1 the lobes have exponen- 
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tially small overlap whence the reduced density matrix becomes 
I 



„SR 



(y,y')\ 



JV»1 



pg \ —\y ± \ — \>\ —\y ± \ — 



(12) 



where pa is given by Eq. 10} . 



III. ATOM-FIELD ENTANGLEMENT 

The RDMs of the atoms, Eq. Q and Eq. are 
derived from the pure ground states in the normal and 
SR phase, respectively. The atom-field entanglement is 
therefore determined by the von Neumann entropy 

S = -Tr(plog 2 p) (13) 

with p the RDM of the atoms (an identical result is ob- 
tained with the field RDM). We present first numerical 
results (already discussed in 01) and a new perturbative 
result for finite N , and then recap our exact solutions in 
the thermodynamic limit. 

We diagonalise the DH in the Fock-Dickc basis, and 
obtain the RDM of the atoms. This is diagonalised and 
the von Neumann entropy is obtained from 

S(p) = -^2 Pk \og 2 p k . (14) 

fe 

where p k are the eigenvalues of the RDM E|. In Fig. ffl 
we plot the results of these numerical calculations. 
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FIG. 1: The von Neumann entropy S between the atomic and 
bosonic modes for a range of system sizes, N = 8, N — 32, 
and N — » oo. The divergence at A/A c = 1 can be clearly 
seen, as well as the strong coupling S = 1 limit. The smaller 
graphs depict the linear entropy Sn n between the atomic and 
bosonic modes Eq. 12611 . and the inverse participation ratio 
P^ 1 Eq. I3UH of the ground state wave-function. The linear 
entropy exhibits the similar behaviour to the von Neumann 
entropy, except the asymptotic limit changes with N because 
of the normalisation. The participation ratio drops to at 
the critical point, indicating a massive derealization of the 
ground state. 



A. Perturbative Results 



/From Rayleigh-Schrodinger perturbation theory, we 
find an iV-indcpcndent result for the von Neumann en- 
tropy for low coupling (with a = X/(u> + u>o)) 

s = -TT^ l0g2 (iT^) - if^ log2 (if^) ' 

(15) 

This matches the numerical data well for A/A c < 0.4. 
Similarly, following Refs. pBll^ for A — > oo, we can iden- 
tify the strong coupling limit ground state as I^gs) = 

75 (l^ -i*> + I - ^>i*>) . wh ere | ± T j x ) is 
a product of a coherent state for the field and an eigen- 
state of J x for the atoms. As this is effectively a maxi- 
mally entangled state of two two-level systems, S — > 1 as 
A — > oo. 



B. Thermodynamic limit 



For N — > oo, the excitation energy diverges as e_ oc 
|A C — A| ziy and the characteristic length diverges as Z_ = 
—j= oc |A — A c | _!/ with the exponents z = 2 and v = 1/4 

|14| . We now use the thermodynamic limit RDMs found 
above to obtain an exact analytical expression for the 
entropy S as N — > oo. This calculation proceeds via 
comparison with the density matrix of a single harmonic 
oscillator of mass m, frequency Q at temperature T [ToL 
|27| with that of our reduced atomic system in Eq.(0}. We 
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find. 



cosh (3Q 



2e_e A 



(e_ — e+) 2 c?s* 



(16) 




2e-e A 



(e_ — e+) 2 c 2 s 2 
(e_-e + ) 2 cV 
2K 2 (e_c 2 + e+s 2 ) 



(17) 



where /3 = l/fceT. We have two equations linking the 
four parameters of the atomic RDM w, wo, A, k (where k 
is the squeezing parameter we introduced in Eq.Q) and 
the three effective parameters of the thermal oscillator 
/?, f2, m. By setting one energy scale of the original sys- 
tem such that wo = 1, and that of the thermal oscillator 
such that m = 1, O = uj, we can uniquely define the 
correspondence between the two systems. 

The squeezing parameter k introduced into the RDM 
in Eq. (0J compensates for the one-mode squeezing that 
the atomic ensemble undergoes as a function of A [l5| , al- 
lowing us to keep the frequency of the thermal oscillator 
constant. With this relation between the parameters of 
the two RDMs, the effective temperature becomes the pa- 
rameter describing the degree of mixing in the RDM. In 
other words, the interaction of the field with the atomic 
ensemble is such that, from the point of the atoms alone, 
it is as if they were at a finite temperature, with the 
temperature given by Eqs lfT7|) . The determination of 
this temperature is not unique, since there are more free 
parameters in Eqs l|17|) than constraints, but the choice 
made here is physically appealing, with the frequency 
of the thermal oscillator constant and the temperature 
varying with A. 

The behaviour of the temperature T with A is shown 
in Fig. |21 and the divergence at the critical point is im- 
mediately obvious. We also plot the squeezing parameter 
k, which vanishes at A c in accordance with the dereal- 
ization of the system here. 

The entropy of a harmonic oscillator at finite temper- 
ature is a standard result from statistical physics 
(setting H = fee = l)j 

S= (^coth(^)-ln(2 S inh(^)))/ln(2). (18) 

Note that this is independent of k, and thus the above dis- 
cussion does not affect the result for S. Solving Eqs l(T7|) 
for the effective parameters, we obtain the von Neumann 
entropy of the atom-field system in the normal phase, 
which is plotted in Fig. We clearly see a divergence 
at A c . 

Moving into the SR phase, if we calculate the entropy 
of a single displaced lobe, exactly the same calculation 
as in the normal phase applies, except with the SR pa- 
rameters instead of the normal phase ones. Around the 
critical point, the entropy diverges and then falls to zero 
for large coupling (not shown here). This is the correct 



scenario for N — > oo, where parity symmetry is broken 
and the system sits in either of the displaced lobes. 

The more interesting case occurs for large but finite N 
where our numerical results indicate that for large cou- 
pling, the entropy S does not tend to zero, but rather to a 
finite value. This can be easily understood by calculating 
the entropy of the positive-parity two-lobcd SR RDM of 
Eq. (|12f> . rather than the broken-parity single- lobe wave 
function as above. We recall that the two-lobe RDM for 
N 3> 1 formally turns out as a mixture of the density ma- 
trices representing each lobe, cf. Eq. (|12f> . A standard 
result [Tl| is that for a density matrix p = \\p\ + A 2 p 2 , 
Ai,2 > 0, the concave nature of the von Neumann entropy 
allows us to write the following inequality, 



S(p)>X 1 S( Pl ) + X 2 S(p 2 ). 



(19) 



The entropy of a mixture of density matrices is also 
bounded from above by 

S(p) < A 1 ^(p 1 ) + A 2 5(p 2 )-A 1 logA 1 -A 2 logA 2 . (20) 

The final two terms are known as the mixing entropy 
[fl| . an d m the case that the ranges of the p\ and p 2 
are pairwise orthogonal, this upper bound becomes an 
equality. Returning to our positive-parity SR RDM, the 
entropy of each of the two lobes is identical S(pi) = S"(p 2 ) 
and they are weighted in an equal superposition Ai = 
A 2 = 1/2. Furthermore, the two lobes are orthogonal, 
and thus from Eqs I|19I2(J|) we have 



S(p) = S( Pl ) + l. 



(21) 



We emphasize that the SR phase entropy plotted in 
Fig-© is a consistent entanglement measure based on 
the underlying pure ground state Eq. 1|11|) . yielding the 
correct large- N behaviour for strong couplings A. 

From Eq. (|18|) we see the entropy depends only on 
the ratio fl/T. In the limit A — > A c , we have Q/T ~ 
y/eZ. Since as A — > A c , e_ — > 0, we see that the effective 
temperature diverges T — > 00 and so does the entropy, 
S — ► 00. As discussed in [TJj, in the neighbourhood of 
the critical point, we have 



1 



1 6A\^.uj i 
"4 ( 16Aj + w 4j ~ ~4 

log 2 I A e — A| + const. 



■ln|A r 



A|]/ln(2) 
(22) 



The prefactor to the logarithmic divergence is identical to 
the exponent characterising the divergence of the length 
scale v — 1/4. Thus we see that, as adjudged by the 
atom-field entropy, the system is critically entangled. 



C. Linear Entropy, Participation Ratio, and the 
Average Linear Entropy 

An alternative measure of entanglement is the linear 
entropy, given by 



L = r,[l-Tr(p 2 )}, 



(23) 
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FIG. 2: Effective entanglement temperature T vs. \/\ c for 
N -» oo, u = w = 1, from Eq. JTIjJ with /3 = 1/T. The 
divergence at the critical point causes the simultaneous di- 
vergence of the von Neumann entropy, signalling 'maximal 
mixing'. The second graph illustrates the behaviour of the 
squeezing parameter k vs A/A c from Eq. 1171 . which tends to 
at the critical point. 



where p is the reduced density matrix of one part of our 
bipartite system, and rj is the normalisation rj = 1+(1/N) 
which gives the correct < L < 1 behaviour [2^ • While 
it is a valid monotonic entanglement measure, it lacks 
some of the full physical i nterp retation provided by the 
von Neumann entropy 0, l29j |. Again, we calculate ex- 
plicit analytical expressions in the thermodynamic limit 
by employing our co-ordinate space ground state (noting 
for D — ► 00 , 77 — > 1 ) , 



Tr(p 2 ) = / dxdx'p(x,x')p(x',x) (24) 

= J dxdx'dydy'ip{x, y)tp(x' , y)ip(x, y')tp(x' , y'). 
In the normal phase 



L 



(e_c 2 + e+s 2 ) 



[(£_ S 2 +6 + C 2 ) 2 (25) 



(e_ S 2 + e + c 2 )(e_-e + ) 2 C 2 S 2 : 



On resonance when to 
simplifies to 

L = 1 



(e_c 2 + e+s 2 ) 

Uq = 1,C 



/2 



(e-+e + Y 



s = 1/V2 this 



(26) 



which is zero at zero coupling, and unity at the critical 
point. In the SR phase we recall the ground state (for 
large but finite N) is a superposition of two lobes, and 
the RDM is a mixture, thus, 



1 



Tr(p 2 ) = -(Tv(pi) + Tr(p 2 ) + 2Tr( Pl p 2 )). 



(27) 



As before, the 
Tr(p 2 ) = Tr(p|), 
we need 



two lobes are pairwise orthogonal, 
and the cross term is zero. Therefore, 



L=l--Tr(p 2 ). 



(28) 



The explicit expression for this is the same as in the nor- 
mal phase, but with the above factor 1/2, and the ap- 
propriate SR parameters. In the large coupling limit, 
Tr(p 2 ) = 1, and thus the linear entropy tends to a con- 
stant 1/2. This function, and the finite numerics, are 
shown in one of the insets of Fig. ^ 



1. Inverse Participation Ratio 

There is also a connection between the linear entropy, 
and the inverse participation ratio 0, a measure of the 
derealization of a wave function. The un-normalised in- 
verse participation ratio is defined as, 



dxdy\^{x,y)\, 



(29) 



Typically, this is normalised over the volume of the co- 
ordinate space, however we work with the un-normalised 
value for convenience. P^ 1 — > for a state delocalizcd 
across the entire co-ordinate space, and P _1 remains fi- 
nite for a localised state, depending on the basis chosen. 

We can interpret the participation ratio as a measure 
of the spread of a wave function over a particular basis, 
akin to the way the entropy is a measure of the spread of 
a density matrix over its diagonal basis. In the normal 
phase, we perform the Gaussian integrals with respect to 
the spin-boson co-ordinates to obtain 



p -i = V^T 

2ir 



(30) 



Thus, in this representation, the participation ratio is 
equal to the Gaussian normalisation factor of the ground 
state, telling us the relative volume in co-ordinate space 
the state occupies. 

In the SR phase, v I / G R |Ar>i = ^(ipi + 1P2), where ^1,2 
represents the two possible displaced lobes. Again, using 
the fact that there is no overlap between these lobes, we 
see, 



P- 



1 f 1 

dxdy-(i[)f + V4) = / dxdy-ipf(x,y), (31) 



thus yet again we can use the normal phase result, with 
the SR phase parameters. This analytical result is plot- 
ted in Fig. where the dclocalization at the critical 
point is clearly shown. 



2. Average Linear Entropy 

To conclude our discussion, we define the average lin- 
ear entropy over all subsystems (iV-atoms and the field 
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In the thermodynamic limit 



FIG. 3: The average linear entropy for system sizes N = 8, 
N — 16, and N — » oo Eq. 1321 which includes contributions 
from both the field and atomic modes. Inset: The derivative 
of the N — * oo limit. 



mode) as, 



Q 



l 



N + 1 



JV-l 

E 



AT 



iV + 1 



1 



AT+1 



(it 



and we have L 



(Jz 



= (dU) 2 

(1 - M 2 ), 



(J,) 2 = (tfb) 2 - N(tfb) + 
= 0. In the super-radiant 
- Nix{dU) + N 2 p 2 /4, and thus 
(fj, = A 2 /A). In both phases 



A<A 
k.N—>oo 



N 2 /4, 

phase, 
r a>a c 

^k,N-» 

the contribution from the mode becomes negligible, and 
Qn—>oo = Lfc- Numerical results for this quantity are 
plotted in Fig. 1(3} . 

It has been shown elsewhere that the average lin- 
ear entropy Q is related to a measure of multipartite 
entanglement proposed by Meyer and Wallach The 
concept of m ultip artite entanglement is a difficult and 
open one (e.g.|33,|3J,|33,|3a|), and while the average lin- 
ear entr opy is limited in the multipartite states it can 
classify j3j it is both easy to calculate and has proved 
useful in a variety of contexts contexts [28l 133 . l37l l38j . 
For completeness, we provide a definition of Meyer and 
Wallach's measure in the appendix. 

In our results Fig.|(2J), we see a clear discontinuity in 
Q between the two phases. This follows directly from 
the discontinuity of the atomic inversion at the critical 
point |l5j , and the simple nature of our symmetric atomic 
states. As a final point, the derivative of Q in the SR 
phase TD limit is -gy = 4A^/A 5 , which we plot as an in- 
set in Fig. (PJ). In the thermodynamic limit the average 
entanglement of all the physical subsystems vanishes in 
the normal phase, but becomes non-zero in the superra- 
diant phase. However, the maximum of Q is not at the 
critical point, contrary to the bipartite partitions and the 
pairwise partitions [lfj. 



where the sum is replaced because the Dicke states are 
symmetric with respect to interchange of atoms. Lk is 
the linear entropy of atom k, and Lb is the linear entropy 
of the mode discussed earlier, 

L k = 7? 2 (1 -Tr(p 2 )), L b = W+1 (l-Tr(p 2 )).(32) 

As before in Ea. (|23*|) . the quantities 772 = 2 and 77^+1 = 
1 + (l/N) provide the correct normalisation. Since the 
linear entropy of one subsystem of a pure state is an en- 
tanglement monotone pHj (it does not increase under lo- 
cal operations and classical communication) , and a valid 
entanglement measure, a concave function of the linear 
entropy is also an entanglement monotone. The average 
linear entropy Q, as we have defined it here, is a concave 
function, and is thus itself an entanglement monotone. 

We can express pk, the reduced density matrix of any 
atom from the ensemble, in terms of the collective expec- 
tation values, 



Pk 



1(1 



(J+i 

N 



,N 



(J-) 

Ifl 4- 2 W ' 
2' 1 T N . 



(33) 



IV. DISCUSSION AND CONCLUSIONS 

Table [3 presents the most important results for the en- 
tanglement measures we have calculated here (and the 
concurrence discussed in |lfj|). and where appropriate 
their derivatives. In particular, we point out the im- 
portance of the divergences at the critical point, and the 
finite size scaling exponents (not discussed here) we cal- 
culated in |lCj . and which have been recently confirmed 
in El. 



fW 


/(A->A C ) 


N Scaling 




\K-X\ 1/2 




l- 


\\c-\\- 1/4 




s 


-Ilog 2 |A-A c | 


log 2 iV(°- 14±0 - 01 ) 


Cr 


1 \/2 
± 2 


jy- 0.25±0.01 


ac R 
a\ 


|A C -A|"i 





TABLE I: Entropy S, Eq. JTHJ, concurrence Cr (from 0), 
and the average linear entropy Eq. 13411 . in the Dicke model 
near the critical point A — > A c . 



and we find Tr(p 2 ) = ± + - 
linear entropy of a single atom is, 



2<J ^i J+> . Thus, the 



L k = 1 - 4(J,) 2 /W 2 



(34) 



For the atom-boson partition, we calculated the en- 
tropy exactly in the thermodynamic limit and numeri- 
cally for finite N. The entropy has a divergence around 
A = A c , which follows the power law divergence of the 



8 



correlation length /_ oc |A — Ad" 1 / 4 . 

There is also a correspondence between the divergence 
of the spin-boson entanglement (entropy) and the der- 
ealization of the wave function. This is highlighted by our 
results for the behaviour of the participation ratio, which 
shows that the ground state of the Dicke model under- 
goes a massive derealization at the critical point. Since 
derealization is a common property of wave functions in 
a quantum chaotic system, our results help strengthen 
the understanding of the relationship between entangle- 
ment and the underlying integrable to chaotic transition 
present in the Dicke Hamiltonian |l5l| . However, most 
generic features of this relationship are still unknown. 
The future of this field lies in closer examination of the 
underlying semi-classical behaviour in quantum systems, 
such as supercritical pitchfork bifurcations ITr| in co- 
ordinate space, or phase space. 

Similarly, we calculated the average linear entropy, 
and example of the Meyer- Wallach multipartite entan- 
glement. Like our other measures, this displays a clear 
discontinuity at the critical point. There is obvious fu- 
ture research to be done in applying the many larger 
classes of multipartite entanglement measures 311 m a 
compact way, and perhaps gaining a clearer understand- 
ing of the behaviour of multipartite entanglement. Other 
avenues of future research may arise from investigating 
quantum phase transitions in other spin-boson models 
[la jig . In particular, while the Dicke model has a nat- 
ural feature that allows infinite system sizes to be inves- 
tigated, there is no reason other spin-boson models with 
non-commuting energy and interaction terms which do 
not have this integrable limit should not exhibit similar 
critical entanglement, with chaotic transitions and level 
statistics. 



system of 7V-qubits, and show its connection to the av- 
erage linear entropy. 

We can write the pure state of N qubits as \ip) = 
6m °6i,— ,b N \bi> •••! friv)- Meyer and Wallach defined 
two new unnormalised states \u k ) and \v k ) as vectors 
in C 2N ~ 1 which are obtained by projecting the original 
state \ip) onto the two possible subspaces spanned by the 
two possible states of the k — th qubit, 



If) = |0> fc ® \u k ) + \l) h ® \v k ). (Al) 
In the Schmidt decomposition, these two subspaces are 
orthogonal (u k \v k ) = 0. Q itself is defined as, 



N 



J2D(\u k ),\v k )) 



(A2) 



fc=i 



where D(\u k ),\v k )) = £ 
alised wedge product 



i<3 \ U i V 3 



u k v k \ is the gener- 



Brcnncn proved that each term D in the sum of Q was 
equal to the linear entropy of the fcth qubit, and thus 
Q(\i/j)) is equivalent to the average linear entropy of all 
the qubits, 



W»=2 



JV-l 



fc=0 



(A3) 
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APPENDIX A: MEYER- WALLACH 
ENTANGLEMENT 

Here we include the definition of the Meyer- Wallach 
entanglement measure |T3 [. as originally intended for a 



where pk is the reduced density matrix of the kth qubit. 
Thus, Q has the required properties of an entanglement 
measure < Q(\ip)) < 1, QQip}) = for product states, 
Q(\%b)) = 1 for the reduced density matrix of every qubit 
being maximally mixed, and Q(\ip)) is invariant under 
local unitaries, both because D is invariant, and Tr(p|) 
is invariant. 

Many pure states fulfil the requirement for maximal 
qubit mixing, and thus give a value of Q = 1. For ex- 
ample, the 3-qubit GHZ state (|000) + |111))/V2 gives 
Q = l, while the Werner state gives Q((|100) + |010) + 
|001))/V3) = §. 
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